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1.0 Quality Compression Optimization 

The following sections deal with the improvement of the comoressed 
mage stated in the report introduction, we bel eve't^ t 

he last significant major problem to be addressed in the s til' 
imagery compression system. 

1-1 Introduction 

The issue is the mechanism which can be used to improve the quality 

cL^r.. ^''''^ ^= parameters^ used in 

compression process. We think this is an important issue, and in our 

If we hold the compression raUo constant, we can compress an imal. 
and obtain on a subjective scale of C. where ^ 

• A is no observable defect, 

• B is observable but not noticeable. 

• C is quite noticeable, but not distracting from the image 

• D IS very noticeable and detracts from the image 

• E IS unacceptable. 

L^'hr°'^"'^''^l '^' P^""'^^"^ by a hand optimization, we find we 
can change the subjective evaluauon from a C to a B with the 

euSantll'r "'^^ "^"if'""'' 'b^t it 

f.H !f ? '^^"'f system which is tailored for each image 

set of parameters irrespective of the image content. 

liZ^Tll °^ compression, consider the following figure. 
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Figure 49 

RGB Image in 24 Bit Color Depth Showing Transform to YIQ 

Figure 49 shows the relationship between the RGB and YIQ color 
represeniauon and the processes , which take place in the algorithm 
as this data is transformed into the final lossy set of wavelet 
coefficients.. The processes which apply to the imagery are shown as 
heavy arrows with the process idenufier shown as part of the arrow 
It no process number is present, then the arrow is just a passive link 
between processes and data. The important point to note is that in 
CO or imagery, we deal with three images (one luminescence, and two 
color planes). The Y component is critical, while the IQ components 
are less sensitive to error introduced by the compression system. 

The processes shown in Figure 49 are as follows: 

• (1) Transform the RGB format into YIQ format, using long 
integers. This format is only an approximation of the YIQ format. 

• (2) Transform the YIQ planes into a wavelet decomposition using 
our own integer wavelet transform. This produces the result 
shown in Figure 50, but for each plane. 

• (3) We have shown an aliernaUve process which is an optional 
down sampling for the IQ color planes. This down sampling may 
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?on^r* '° P'"*^"" ''"^5'= P'^n« cither one 

fourth or one sixteenth the size of the original plane. If this 

process is to be done, it will be accomplished prior to the wavelet 
transform of process (2). r y k> uic wavciet 

(4) Here the first step in the loss occurs. The wavelet coefficients 
are now quantized to a number of levels depending upon which 
quadrant IS being processed, and the desired compression or 
quality factor. 



• (5) Simultaneously with step (4), the wavelet coefficients are 
being matched against threshold values, and if the values are less 
than the established threshold values specified, then the resultant 
value IS set to zero. 

• (6) The last step in the process is to entropy compress the 
resultant coefficients using cither Arithmetic. Run Length or 
Huffman, of Huffman and Run Length combined. The key issue is 
the amount of compresiiion desired against the invested time 
required to get that level of compression. 

The issue now. irrespecUvc of the down sampling or not of the IQ 
components, is the fact that we process the three planes into five 
levels m a decomposition as shown in Figure 50. From this figure 
one can see a total of 16 regions (quadrants) which arc defined by 
the numbers 1 through 16. Each of these sixteen quadrants have two 
additional parameters associated with them. The parameters define 

quantization and threshold values for that particular quadrant 
Since there arc three planes for color (only one for gray level) the 
maximum number of parameters that the u.ser can control is 96 -- 16 
for quantizaUon and 16 for thresholding for each of the three color 
layers. In the case of a gray level image, there are only parameters 
for one layer. 

Our experience has shown the parameters for an image are very 
scnsiuve in some cases, and not in other cases. In order to measure 
this sensitivity, we generated the variance of the wavelet coefficients 
in the 16 quadrants. Table 1 provides these values for each of the 
three planes. 
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Figure 50 

Decomposition of a Plane into Five Levels 
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36 
23 
10 
128 
35 
37 
233 
73 
64 
499 
138 
156 



4.7 
3.2 
1.3 
5.8 
4.3 
2.8 
11.4 
7.4 
3.3 
27.5 
11.8 
6.3 
41.7 
11.2 
12.2 
'2818 



1.4 
1.1 
0.6 
1.7 
I.l 
0.8 
2.6 
1.1 
1.0 
4.1 
1.5 
1.1 
10.1 
2.3 
2.0 
.221L 



164.2 
134.6 
28.9 
276.1 
274.3 
124.8 
223.7 
298.1 
1 14.2 
174.8 
285.1 
107.9 
135.4 
245.0 
89.3 
7-408 



L 

27.3 
24.9 
4.6 
37.7 
67.3 
22.7 
73.2 
107.3 
35.5 
73.9 
128.2 
47.8 
75.6 
144.0 
49.1 
6904 



6.6 

3.0 
1.1 
12.9 
12.4 

5.0 
13.7 
16.2 

6.B 
10.7 
14.9 

6.7 

8.7 
11.2 

3.1 

77,9 ,.„ 



Table 1 

Wavelet Coefficient Variance by Numbered Quadrant for Two Images 
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The clear information to be gained by the numbers in Table I is the 
images are quite different, and the quantization or threshold levels 
set for all images will only be an approximate solution at best. The 
optimum solution would be to have the quantization and threshold 
values set according to the variance values. Such settings could be 
found in a table with various ranges, and for each such range, the 
parameters of interest could be defined. This would give a more 
optimal, solution, but still not the optimal solution. In order to get 
optimality, one would need to search over the variable space using 
the near optimal settings to find the actual best values for the 
parameters. 

As the values in Table 1 grow, the implication is a finer mesh sire for 
quantization. That is, the number of bits one wishes to allocate to the 
output could be varied by quadrant. Those quadrants with large 
variances will utilize more bits, while these with low variants will 
uUlize fewer bits. In this way, the number of bits resulung from 
quantization will remain the same, but their allocation will differ 
depending upon the nature of the image. 

1.2 Approach 

The solution to the problem posed in the previous section is to 
determine how the ninety-six parameters interact with one another. 
The problem is a bit more sophisticated than just measuring 
parameters, however. The issue is with each parameter change, the 
compression ratio will change. If a compression ratio, or a quality 
factor which indirectly defines a compression ratio, is specified, then 
the user wants the compression ratio to remain identical over the 
changes in the parameters. In order to accomplish this, there are 
two parameters which we must monitor: PSNR (peak signal to noise 
ratio which is defined to be PSNR = 20 logio (X/MSE) where the X is 
the average absolute value of the pixels in the after image and MSE 
is the mean squared error measured between the before and after 
image) and the compression ratio. The compression ratio must be 
held constant, and the PSNR needs to increase, and the way to 
increase the PSNR is to reduce the MSE. 

The difficulty with this system as described is in many cases, small 
changes in the parameters introduce significant changes in the MSE. 
Also, we believe, the parameters are not independeot. We have also 
seen images where the parameters can be changed in one way, then 
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altered, and the results are exactly the same. This indicates the 
opUmal value is not a single point, but rather something like a plane 
with little slope. " *^ 



1.3 Status 

We have established the rules under which the optimal solution will 
need to exist, and arc at the moment writing software to measure the 
variance within the Lightning Strike environment. Once this is done, 
we will being examining many images to sec how close we can 
determine the optimal parameters for a defined set of variances. 



• 
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CIS-2 Image Compression Algorithm 



HONGYANG CHAO 



Part I: Brief Review of LSIC 3.0 
L Introduction 

CTS-2 (tcmporaiy name), which has been being used in Lightning Strike 3.0 image 
compression software, is a wavelet based image compression algorithm. CIS-2 has 
following inventions: 

• Integer reversible wavelet algorithm with Property of Precision Preservation; 

• Subband oriented quantization and related entropy coding; 

• Wavelet lossless compression for color and gray images; 

• Progressive transmission algorithm for color bit compression; 

• Progressive transmission and decompression algorithms; 

• Non-uniform image compression algorithm; 

• Quality based wavelet coefllcicnt quantization tables; 

• Attached optional post-processing filters; 

• Image map editor; 

• Optional peak signal noise ratio controlled comprcvssion; 

• Special split and merge wavelet compression algorithm for very big image 
compression without any boundary efifccts ; 

• Image dependent parameter optimization • 

2. Main steps of the algorithm 

In Lsic 3.0, three kbds of different image compression methods arc included: 



Section 2. 1-2 J will give brief description of above method The details will discuss 
later. 

2.L Main step of method 1 

Figure 1 and Figure 2 give the flow charts of the image compression and 
decompression of method! respectively. Every step in both compression and 
decompression has lot of details, which will be described later. 



Method 1 : Quality controlled wavelet based compression 
Method 2: Color bit depth compression 
Method 3: Wavelet lossless compression 
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Figure 1: Compression flow chart for Method 1 
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Figare 2: Decompression flow chart for Method 1 
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2.2. Main step of method 2 

This method based on using less number of colors to approximately represent 
original images. Following figure gives the main step of the algorithm for the 
compression and decompression. 



Input Original Image 
Accept the Number of 
Color 




Create Color Table by Color 
Quantization or Input a 
Fixed color Table 





Wavelet Transform 
on the Index 




Calculate the Index for 
Every Pixel with Dithering 










Entropy Coding the 
Transformed Index 




Output the Coded File 





Figure 3: Compression flow chart for Method 2 
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Fig:ure 4: Decompression flow chart for Method 2 



23, Main itcps of method 3 

The main steps of method 3 is almost as same as method 1. However, at the every 
step we use different methods. Following are its compression and decompression 



flow charts: 
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Figure 5: Compression flow chart for Method 3 
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Figure 6: Decompression flow chart for Method 3 



3. Brief Description for Other Features 

Most Features (or inventions) are included in above three compression methods. 
Following is the brief introduction for some inventions list above. 
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3.1, Progressive decompressioa 

This aJgorithm allows users to decode images from the lowest resolution to 
highest resolution. The advantage of this feature is that usen can dosvnload small 
piece of the coded file and view the image at lower resoluUon to determine if they 
want to download the whole image. 

Steps: 

(a) Input the lowest pass component U' of the coded file and reconstruct the 
lowest resolution image /' ; 

(b) Display image /'.If the user doesn't like it or the resolution is big enough for 
, stop; otherwise, go to next step; 

(c) Input the lowest three band-pass components //I", iZ/'ahd/f//" 
successively in the current coded file. Reconstruct the new image /' from 
LL' . HL\ LH' and HH° . Let 7" - /' , go to step (b). 

3^. Non-unifonn image compression 

The algorithm allows users to divide the image into several parts with different 
interests and compress these areas with different qualities. The areas can have any 
shape. This algorithm is only available for method 1. 

Original image goes Lhough all of the procedure except quantization part, which 
follows the steps below: 

(a) Creating the bitmap matrices related to the areas chosen by the user; 

(b) Wavelet transform to every bitmap matrix; 

(c) Different quantization in different areas according to the transformed matrices 
obtained above step. 

3 J. Peak Signal Noise Ratio (PSNR) controUcd compreuion 

Peak Signal Noise RaUo (PSNR) is an image quality measurement used by most 
professional people. PSl^ controUed compression aUows users to choose their 
desired PSNR for the compressed image. 

The related algorithm is an iterated system: 

(a) Picking an initial parameter setting /J ; 

(b) Quantize the wavelet coefficients with and calculate the corresponding 
PSNR; 

(c) If the PSNR is close to desired one, stop and output the coded file; otherwise, 
get an adjusted vector and set ?, ^ ?, + A?, , go to step (b); 
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3.4. Attached optional post-processing filten 

Users can choose any number of follovnng processing filters at their compressing 
time. The desired results can stored in the coded image file, and. anyone who 
decompress the coded file will see the same result immediately. 

• Sharpening images 

• Smoothing images 

• Improving the visual quality 

• Brightening the images 

Image map editor 

Image Map Editor create a image map over Lsic.?.0 compressed image file. User 
selects one or several areas of compressed image, assigns the http links to the areas, 
then. Image Map Editor calculates the coordinates of the areas and outputs a HTML 
associate with the image. User can add such information into hisAier source code. 

Following is an example of such image map: 
C^<EM3ED SR6^cow.C43d" type = "image/cis<od" WIDTH- "237" poiy= "44, 
45, 103. 78, 103, 86, 5^86, 54, 78", href= •'http://www.infinop.com"X/EMBED> ' 

3.6. SpUt and merge wavelet algorithm for very big image conpressioa 

This algorithm allows users to compress very big image by an ordinary machine. 
The key is to divide the original image into several smaller pieces and 
compress/decompress them separately by using overlap and dis-overlap technique. 
With this technique, the compression/decompression piece by piece is equal to 
compress/decompress the whole image together, which means users won't see any 
edge effect at decorapirssed image which appears at general split method. 

Also, with this algorithm, users can cither decompress the whole image or choose 
the specific part to decompress according to a image map we create for the division . 

3.7. loaf t depcsdeat optimized param eter setting 

This algorithm allows user to gel the best (or almost best ) image quality at the 
desired compression ratio by choosing image related parameter setting. 

3 J. Integer reversible wavelet algorithm with PPF property 

See attached unpublished paper tiUed as "An Approach to Fast Integer Reversible 
Wavelet Transformations for Image Compression" 
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3.9. iDtrgtr color trauformstioD 

This algorithm is a integer reversible transform which has been used in lossless 
color image compression (Method 3) for Lsic 3.0. 

The algorithm transform RGB color components to a new set of color 
components Y-Nb-Nr: 

Forward transform RGB to Y-Nb-Nr: 

r=c;+int{Y). 

A*-5"Int(f), 

A^rBif-int(f). 
Inverse transform Y-Nb-Nr to RGB: 

G = r-Int(^). 

3.10. Sabband rebittd QaaotuatioQ mad entropy coding 

rnis entropy coding method is just designed for wavelet based image 
compression. The main idea is to take the advantage of difFeient quantization at 
different subbands and encode each band according to its content This method 
reduce the coding cost greatly. 



CIS-1 Image Compression Algorithm 



HONGYANG CHAO* 
Computer and Information Science Inc. 

1. Introduction 

CIS-1, which has been being used in Lightning Strike image compression software, is a 
wavelet based image compression algorithm. CIS-1 has following advantages: 
0 Reach almost optimal compression ratio; 

0 Keep the major characteristics as more as possible. In other words, it reduce 

insignificant components gradually according to human visual system, so that people 
can still accept the image quality at the extremely high compression ratio; 

0 Fast 



2. Main steps of the algorithm 

Figure 1 and Figure 2 give the flow charts of the image compression and decompression 
respectively. Every step in both compression and decompression has lot of details, which 
will be described later. 
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Figure 1: Image compression 
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Figure 2: Image decompression 

3. Image Compression 

For the compression, we only have to describe three parts: image transformations 
quantization and thresholding, and entropy coding. 

3.1. Image Transformatiou 

TTie image transformations involved in this algorithm include color transform (for color 
images) and wavelet decomposition ( for both gray level images and color images). 

(a) Color transform 

In general, input color images are based on RGB color model, such as 

TIFF or BMP images. In order to get high compression ratio, it is better to change RGB 

color model to other color models, such as YIQ or YUV models 

RGB to YIQ: 

[rl ro.299 0.587 0.114 Ti?l 
I / 1 = 1 0.596 -0.275 -0.321 1 G I. 
LfiJ Lo.212 -0.523 0.311 i^J 

RGB to YUV: 



frl ro.299 0.587 0.1 Mji?! 
I C/1 = 1-0.148 -0.289 0.439 I G I. 
W Lo.615 -0.515 -0.1 i^J 



(b) Wavelet decomposition 

The purpose of wavelet transform is to represent the original image by different basis to 
achieve the objective of decorrelation. There are a lot of wavelets which can be used in 
this step. In CIS-1, we use a wavelet which results in the following algorithm: Suppose 
C**=[c2t]^j^ (J =0,-\M-U * = 0,-".//-l) isoriginalimage. whereA/and// are 
integers which have the common factor 2^ (I is a positive integer). After one-level 
wavelet decomposition, we will get four parts as shown in figure 3. 




Wavelet 



Transform 





HD^ 


VD^ 


DD^ 



Figure 3. Wavelet image decomposition 



We caU = [c]j] (;• = 0, " s - 1 ; * = 0,. . - 1) the blurred image of , HD^ the 
horizontal high frequency part of VD' the vertical high frequency part, and DD* the 
diagonal ones. Setting -C\ we can repeat the same procedure L times or until the 
size of the new blurred image C' is small enougL Therefore, we only have to give the 
algorithm for one level decomposition: 

(1) Let C** - rC** , r >0 is a factor which can be changed for different needs. 

(2) Transform for image columns: 
0 For *aO,---,^-l, calculate 



"o* ~ 2 ' 
1,- 



(3.1.1) 



o For * = 0,* 1 , calculate 



(3.1.2) 



(3) Transform for rows: 

0 For y = 0, • . <^ - 1 , computing 



and 



(3.1.3) 



I ^JO - 



3 '* 1 



I 

0 For = 0,- • - 1 , computing 



and 



(3.1.5) 



j<-^.u.-^^^.* = l.-...^.2. (3.1.6) 
(4) C - [cj4 HD' » VD' = [v^,,] and DD' = [daj J. y = O.-,-- - 1 

* = o..-.f-i. 



Remart. If it is necessaiy, we also can use matrix multiply 

Wavalet CoefBcient Image of / levels = W/C°Wj 



Here, is the transfonn matrix for / level wavelet decomposition. 
3J. Threstaoldiiig and Quantization 

Both thresholding and Quantization aUow us to reduce accuracy with which the wavelet 
coefficients are represented when converting the wavelet decomposition to an integer 
representation. TTiis can be very important in image compression, as it tends to make 
many coefficients zeros-especiaUy those for high spatial frequencies. 

After L level, for example 1=3. ^vavelet decomposition, we get the wavelet coefficients 
of the original image as plotted in Figure 4: 





HD^ 
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HD ^ 
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DD ^ 



Figure 4. L=3 wavelet coefficients distribution 

(a) Thresholding 

In algorithm CIS-1. we use multilevel uniform thresholding method: Let 

be the chosen thresholds, where /, is the threshold for /th (/ = \,.-.L) level and / is a 
threshold for blurred image . Hiresholding is to set every entry in the blocks <?', 



HD', VD' and DD' (/ = I, to be zeros if its absolute value is not greater than the 
corresponding threshold. 



Remark For color image, we can have three threshold vectors which correspond three 
different color bands, such as Y, I and Q. 

(b) Quantization 

Quantization is to scale the wavelet coefficients and truncate them to integer values. In 
CIS-1, we use the quantization table shown in Table 1 to implement it 



_I 


_2 


• • • 






„l 
IVD 


„1 


• • • 






Qua 




• • • 


L 

<fnn 





Table 1. Quantization table 
Here, the entries ?^ are quantization factors for blocks HD' (/ = 1, • • , q^yo and 
for blocks and DD' = -..I) respectively, and ifae factor q^*' isforthemost 
blurred image . All the factors are integers between 0 and 255. The quantization 
scheme for the block HD' (/ = I, is 

U*=round{ ^?'° ).y°0.-..^-l; ;t = 0,- -.f - 1. (3.2.1) 



HD 



Here, hJj^ = 0, • • ^ - 1; * = 0,- - - 1) are quantized wavelet coefficients in block 
HD' = .,1), 



maxx/n = max 



and the function round (x) gives the nearest integer of x . The scheme of quantization for 
other blocks are the similar to (3.2.1). 

Remark. For color image, as the same as thresholding, we can have three separate 
quantization tables for different color bands. 



33. Entropy Coding 

Here, the encoding means the lossless compression for the wavelet coefficients. It is 
divided into two parts: Coefficient arrangement and entropy coding (Hufi6nan or 
arithmetic). 



4. Decompression 
4.1 Decoding 

Decoding, just as encoding, can be divided into two parts: Entropy decoding (Huf&nan or 
arithmetic), and coefficient rearranging. 

i2 Dequantizatioo 

After Decoding, we get quantized wavelet coefficients in i ♦L+i Blocks. Dequantizing 
uses the same quantization table as quantizing, and the scheme as follow: for / = 1,- • 

(4.2.1) allows us to get the approximate coefficients for the blocks HD' (/ = U-^Z) , 
which is shown in Figure 4 . The dequantizing scheme for other blocks are similar to 

(4.1.2) . 

4 J Inverse Image Transformations 

(a) Wavelet reconstruction 

We are going to describe the algorithm for one-level reconstruction which is plotted in 
Figure 5. 



'^•^-^.;=0.-.if-l; * = 0.-.f-l. (4.2.1) 
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Figure 5. One-level wavelet reconstruction 



(1) Inverse transform for rows: 
0 For yaO,--,^-l, calculate 



and 



L, , / 

I 



(4.3.1) 



0 For ;=0,"-,^-l, calculate 



1-1 1 



I 



aad 



(2) Inverse transform for column: 

0 For ^ = 0,- • A'' - 1 , calculate 



1 



|4 = 4+ 2 ' 



and 



OA' 



(4.3.2) 



(4.3.3) 



(4.3.4) 



(4.3.5) 



(4.3.6) 



(3)cJ^=c;^/r. y=0,...,3/-l; * = O.-.iV-l. C" =[c;,]^^. 
(b) Inverse color traosfonn 

For color image, we have to do inverse color transform 
0 YIQ to RGB 

[R] fl-OOO 0.956 0.621 in 
IG 1 = 1 1.000 -0.272 -0.64711 / I. 
UJ Ll-OOO -1.106 1.703 JLeJ 



o For y = 0,---,^-l, calculate 

|2;..-'^yo+ 2 



1-1 1 



(4.3.2) 



(4.3.3) 



and 



(2) Inverse transform for column: 

0 For ^ = 0,- 7/ - 1 , calculate 

[-0 .1 4)t 

|<^u =^04 + — ; — I 



and 



, 0 



-0 _ 



(4.3.4) 



(4.3.5) 



9 , 20 



(4.3.6) 



y=0.-.A/-l; *=0.-.-.iV-l.C''=[c;_,]^^. 
(b) Inverse color transform 

For color image, we have to do inverse color transform 
o YIQ to RGB 

fl.OOO 0.956 0.621 1X1 
I G I = 1 1. 000 -0272 -0.647 II / I . 
W Ll.OOO -1.106 1.703 JLfiJ 



0 YUVtoRGB 



4.4 Neceuary Pott Image Froceuing 
(a) Color Quadzadon 



fl.OOO 0.000 L140 Ti'l 
1 1.000 -0.395 -0.581 1 C/i 
Ll.000 2.032 0.000 IkJ 
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Abstract 



In this paper, wc propose a general method for creating integer wavelet transforms which can be 
used in both lossless (reversible) and lossy compression of signals and images with arbitrary sire. 
This method allows us to get a scries of transformations which are very close to the corresponding 
biorthogonal wavelet transforms or some non-orthogonal wavelet transforms, but can be cauculated 
with only integer addition and bit-shift operations. In addition, the integer wavelet transforms 
created in this paper possess a property of precision preservation (PPP). This propertv is very 
useful for conserving memory in both compression and decompression, and s-peed up the whole 
procedure in some applicauons. The motivation of this paper comes from the lifting scheme [11 
and S+P transform [3j. 
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1. Introductioa 

The wavelet transform has proven to be one of the most powerful tools in the field of image 
compression. Theoretically, the wavelet transformation is lossi:t>s, but since all computers have 
only finite precision, most of transformations are lossy in practice, even when wc use floating 
point calculations. On the other hand, integer calculations are much faster than floating point for 
virtually all computers; and integer computations are much easier to implement in hardware which 
is more important in some ajqjlications. The memory utilization of integers is also a positive 
consideration. The difficulty is, if we directly use integers in the wavelet transform and its inverse 
without some proper considerations, it will cause the loss of accuracy. For some important image 
applications, the user wants to have complete control of the precision in which the image pixels are 
represented during the compression process, and thus prefers to have the image compressed from 
lossless to lossy. 

Lossless compression is also very important for images found in such applications as 
medical and space science. In such situations, the designer of the compression algonthm must be 
very careful to avoid discarding any information that may be required or even useful at some later 
point. From the academic point of view, it is also very interesting to have a compression scheme 
which has very fast performajK:e, and which can exactly reconstruct the image when necessary. In 
addition, it is also very useful to have some wavelet transforms which exhibit the property of 
precision preservation (PPP), which can be utilized in the computer which has limited precision 
and limited memory without losing any precision during the computation. 

In this p^)cr, we arc going to describe two general methods from which one can get the 
integer wavelet transform desired. All of the wavelet transforms from the methods given in this 
paper possess the property of precision preservation (PPP). We draw on the work of several other 
authors who have already contributed to this area [2-3], where some specific examples were 
developed. However, this paper presents a more general method which allows one to sec several 
new results as well as those presented and acknowledged prior to this work. 

This paper is organized as follows: Section 2 and 3 give some examples of integer wavelet 
transforms. The examples in section 2 arc the starting point for our approach, and the examples in 
Section 3 show the steps and motivation of our general method. Section 4 indicates how one can 
use the lifting technique to create an integer biorthogonal wavelet transform. The Correction 
technique to generate more general integer wavelet transforms is described in Section 5. Section 6 
describes how to process boundaries in order to apply the integer calculation in finite sized irnages 
or signals. In Section 7, we prove the integer wavelet transforms developed by both the lifting and 
correction method possess the property of precision preservation (PPP). Some example images 
are also shown in this scaion. The last section. Section 8, provides the conclusion to this paper. 
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2. Basic integer wavelet transformations 

We provide the following two examples as the j;tarting point for our new method. For the 
sake of convenience, length, and simplicity, we only discuss the algorithm for a one level 
decomposition and reconstruction and only for a one dimensional signal. The extension to two 
dimensions is immediate as the rows and columns can be treated into a sequence of one 
dimensional signals. For the following examples, assume that {c°}2^^ is the original signal where 
the superscript indicates level and the subscript indicates a particular point in the signal. Also, 
{ci and ' arc its decomposition parts at the fiisi level. Heir 

an even number, 
> an odd number ; 

{c\ and }^"' are its low frequency (0 part and high frequency {/?) part, respectively. For 
multi-levels, we just treat {ci ]^^* as and repeat the procedure again. 



ff. ifiVisa 



E««mpU 1: A (2,2)*wavele( transform by integer calculation. 

This transformation is similar to a variation of the Haar wavelet transform which uses low 
and high pass analysis (decomposition ) filters given as: 



(1) Compute 

(2) Compute 



n 


0 


1 




1/2 


1/2 


s. 


1/2 


-1/2 



c; =:Im(^)+cL^, Jfc = 0, -.A/,-2, 

is an even number. 



cj^. if A/ is an odd number. 

Here, Int(jt) is an arbitrary rounding funaion which may have different interpretations. 
For example, Int(x) can be the integer which is nearest to x , or Int(x) may be any integer 

which satisfies x - 1 < Int(x) S z , etc. It is easy to see that all entries in both {ci and 
{^'-Ko * are integers. 



(2.1) 



(2.2) 
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From (2.1H2.2). we can easily get U.e following i„,,g„ ^nsmiction algorithm: 

(b) Reconstruction 

(l)If \ « an even number, compute 

or. if iv is an odd number, we have 

' T X 



<^t.i = - InK-^) . /t = 0.- • •. A', - 2. 



(2) Compute 



(2.5) 



fi^ ^•"-C2.lH2.6)arenot]inearbecauseoftheroundingoperationInt(x) this 

was apphed fim to the coiunuis and then to the rows, the inverse t^,fonn3tioc^t«t be 
appued first to the lows and then to the columns. 

£um{tlL2: Lazy wavelet transform. 

'2^!^^^"' i"-"'^ - important 

'^""P^"^^ « nothing else but sub-sampling the even 

j^^do^ mdexed samples. Decomposition and ™cu^^ 

• Examples 1 and 2 are not good transforms for image compr«sion. but they are simple 

laslirrr'^"'^^^'''"'"^"^^^^- ^«above,weLin 
only as astaituig pouu for our integer. «venible. wavelet transform algorithm 

n,u« mention that there is another interesting prt^erty in the above two transfonns 

winc^maynot^easilyseen.Iftl.v.,uesofthesigna,pixelsa.rep.es«^ 

&^ the rcconnrucuon algorithm, the computer Will get 
trough Che s^^complementary code propeny. We call this property aPraperry^r^ 
i»rMerve«wnfW) for these wavelew. j f^r^ of rrecuion 
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It is known that the general values of the high frequency wavelet coefficients am small, and 
all higher levels in the decomposition also provide generally small values in the high frequency 
band. This allows the preservation of precision during the computational stage of the wavelet 
coefficicnu. Now, the complcmcntaiy code property, the other aspect of the PPP property is .a 
well known characteristic of integer arithinctic as done by the computer. Consider the computation 
of the diffeirncc of two integers given as c = i-a and the inverse computation of a ■ 6 - c. The 
nature of the computation within the computer can be specified as follows: 

b-a if -t^^^^b-aKl^'^-X 

l^-^b-a if 6-fl<-2^-' 



and the inverse is 



a- = 



-2^^b^c„ 



if b^c^^l^^' 
if &-c,<-2^-* 

where the m subscript indicates the internal rspresentation, and the range of the integers a, b. c is 
[-2*"^ 2**' - 1]. The internal representation of when it is outside the range, its appearance is as 
a two's complement nunabcr, so the reprcscnution may not be the same as die external 
representation of c. However, the .«;amc coraplcmcntaiy code for the will cause the internal 
representation to be identical to the externa] representation of a. For example, if we let b^2 
(OOOOOOlO) and a=-/27 (10000001) then c^has the internal binary value of (lOOOOOOl) when 
q^. With a value of -127 for Cm, the inverse value for a„, will just be a. 

In fact, for Example I, this property is obviously true. While for Example 1 , if the range 
of the pixel values is within a finite number of bits, say q, wc can only use q bits as the working 
unit, which means the value of transform coefficients will also be within the interval with length 
2*^, say [-2^', 2^* - 1]. Due to the nature of computation on a machine, most machines will 
implement (2. 1 H2-2) automatically as follows (the complementary code property): 



if c° -cLhS2'-. 



if 



(2.6) 
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lD»(Y) + C:i->. if •2'-'iIm(|-) + c,V,<2'-'. 
Int(y)+4*,-2'. if Tnt(^)+c«^2 2'-'. 

+ if In((^) + cL<-2'-'. 



(2.7) 



While the rcconsmiction algorithn, (2.3) and (2.3) will be implemenied by the computer itself 

.f -r-'^c\-int{^)<r. 



as 



c;? ^ 



ci-lnt(f), 
2'*[c;-Int(|-)^ 

j1 



if <:;-Int(-:A-)<-2»-^ 



c{ - liit(-^) ]- 2* . if d- Int(i) >2^. 



2 

/I j.>.0 



if -2'-'Srf| + 4.,<2'-' 
f« '1<^;+c»,^+2'. if <f; + 4„<-2»-^. 



[<'^fL-2*. if c/;+c;..S2*-'. 



(2.8) 



(2.9) 



It is obvious that (2.8).(2.9) arc just the .tvcr>c of (2.6>(2.7). It i. also easy to see that if wc 
properly take advantage of the bound in the coefficient size mentioned above, the algorithm can be 
implemented using a minimal amount of storage. 

3. More Examples and AddlUonal Analysis 

In this section we are going to give more examples which will give some motivation for our new 
approach. 

Euffi&kJ: A a6)-wavelet transform by integer calculation [2], 



This transfo rmadon is similar to usin^ fol 



0 

-1/16 



owing analysis filters 



•1 



0 

■1/16 



1/2 
1/2 



1 



1/2 



0 

1/16 



0 

1/16 



(a) Decomposition 

Decomposition suutx with Example 1 at step (1) and (2). and then upgrades the high 
frequency component at step (3): 
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(1) Compute 

(2) Compute 



c;»Int(:^) + c«,,. i.o,.,M-2, 



Ini(-^) + e9_^ , if is an even number, 



(3) Compute 



ifA'isanodd number. 



and then, if ^ is even, calculate 



the, calculate 



(b) Reconstnictian 



Tie rtconstruction algorithm is identical to the decomposition algorithm, except 
ninmng "backwards". ^ y 

(I) Compute 



It is now 



rfi;«=Ini(lzl),rf. 

4 



aadtheiuif ^ isevcD, calculate 



else, calculate 



(2) If ^ is an even number, compute 



0 I r dt" 

c,u,=C4'Int(-i-). *sO,...Ar,-i ; 
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or. if is an odd number, wc have 



(3) Compute 



0 I 

Cn*,=c,-Int(-^). i = 0.. ..7V,-2. 



Wc see in step (2>.(3) above, that they are just the same as shown fnr ,1. 

the (2.2>wavelet transform (Exan^plc,). ''''^^ '«°n«n.ction of 

EaoBiU: A (1.3).wavclct cransform by integer calculation , 



(a) Decomposition 



n 


•1 


0 


1 




1 


0 


0 


f. 


1/4 


•1/2 


1/4 • 



frequency coMponent ai step (2): 



(DSet 



= ci. 4=0,-.A/,_l; 



(2) If A/ is an even number, calculate 



j;=Iat(-^).rf« ifc = o...,i»f,-2. 

Otbenwe, if ^ is an odd number. caJcuiate 

^;«Int(fi^).,.^^, i = 0...,^,-,. 

(b) Recoiutniction 

(2) UN is an even number, calculate 

c^. = Int(iU-^)_rf;. * = o...,Af,-2. 
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Othenvise, if // is an cxld number, calculate 



E' l mwlf A (5.3)-wavclct transform by integer calculauon. 

This transfomution is also simiJar in function to using the biorthogonaj analysis filters. It 
is given by 



n 


.2 


-I 


0 


1 


2 




• 1/8 


1/4 


3/4 


1/4 


-l/g 




1/4 


-1/2 


1/4 


0 


0 



(a) Decomposition 

This decomposition starts with Example 3 at step (1) and upgrade low frequency 



components at step (2): 
(1) Set 

If is an even number, calculate 



r 



Otherwise, if is an odd number, calculate 



(2) If N is an even number, compute 



ci = ci»-lnt(l). 



<^'Ni-.=ci».,-Int( 



Otherwise, if /V is an odd number, calculate 
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(b) Reconstruction 
(i) Compute 



4=ci''-Int(^). 

c;=c'.'-Int(^k_t£L). i-l.....Ar .2. 



c«=c;+Int(^). 



Then, if N is even, calculate 



else calculate 
(2) Compute 



-0 I . T,.^ ^^-' x 



Then, if is even, calculate ci_. = , - 

The PPP property for Example 1 -2 mentioned at the end of the previous section is also 
applicable for these three examples. It is obvious these three transformaUons are not reaUy linear, 
but they are simUar to the one using the corresponding filters given above. Especially, the filten in 
Example 3 and Example 5 belong to, with minor modification, the group of the best biorthogonal 
filters for image compression according to both our experience and the conclusion of (4]. 

Also, from the above three examples, we can note that if we begin with integer (linear or 
nonlinear) wavelet transformations and then use some proper upgrading formulas, wc can get 
other, much better integer, wavelet transfonnations for image compression. Now. the key 
problem is: What kind of deductive formulas should be used? We provide an answer to this 
question in the following two sections. Section 4 and Section 5. 

4. Lifting Scheme and Integer Biorthogonal Filtering 

The Ufting scheme, discovered by Sweldens ( I], is a new approach for constructing 
biorthogonal wavelets with compact support. However, the most interesUng part of this method 
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for us is: it can be used, with minor modification, 10 create integer biorihogonal wavelet 
transfonnations. The foUowing is an adaptation of the technique of [1]. 

Definition 1 . The set of filters {h, h, g,g] is ^ set of biorihogonal filters if the following formulz 
is satisfied: 

and Ko»^^h,e'^ and g{€o)^^g^e' 



where m(a)) = 



i 4 



and similarly for m(Q)), h{co) and g{co) . 

The following lemnia is the main result of the lifting scheme [I] reported as corollary 6 in 
that paper. 

Lemma 1 Take an initial set of finite biorthogonal filters [h, A g°, g], then a new set of 
finite biortbogonal filters {h, h, g, g] cdnbc found as 

g{0))=g'i(0)'h{CD)s{2C0)- 

Similarly, if we take {h^,h, g, as an initial sciof biorthogonal filters, anew set of finite 
biorthpgonal filters {h, h, g, g] can be found as 

h{Q})^h\0))^ giQ))Ja^ 

- • (4.1b) 
g{(o)^g\Q))-h{co)7{2(o), 

Here, s{Q)) is a trigonometiic polynomial and the corresponding filter s is finite, and so is 5(iU). 
Actually^ regarding the filters, (4.1) is equivalent to 

0 V. ^^-^^ 

^ ^ ' . (4.2b) 
t 

Next, we use the lifting scheme with minor modification to create an integer, nonlinear, 
quasi-biorthogonal, wavelet algorithm. Suppose {c°} is a original signal, {c^ } and ^\ }arc 
again its low and high frequency decomposition parts, obtained by using the filters [h, h, g, g). 
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If we use fdters {h.g] for decomposition (analysis), the corresponding decompositioa algorithm 



IS 



While the rcconstniction algorithm will be 

^0 



relatc<J to the synthesis filter {h,g\. Here, parameters a. and are positive constant with 
cr, • = 2. For example, in the situation of regulai bionhogonal decomposition and 
reconstruction, = = ^2 ; and for Example 1 through Example 5 above, a, = 1 and = 2 . 

If the set of filters [h, h.g.g] is from {h, k\ g", g } by (4.2b). then the decomposiUon 
can be accomplished as follows: 



Calculate 



2. Calculate 
The relative reconstruction scheme will be 

V 

1 . Calciilare 



(4.4) 



(4.3) 



2. Calculate 



/'.i.o 



0 > 



(4.6) 



Here, equations (4.3) and (4.6) are just the wavelet ( inverse ) transforms using biorthogonal 
filters [K h", g^.g]. While (4.4) and (4.5) are forward and backward upgrading formulas. 

Similarly, iftho set of niters [h, h, g. g] i.s from the initial set of filters [h\ h, g, g"] 
by using (4.2b), the relative decomposition is: 



Calculate 
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2. Calculate d[ = d' " - r 

The reconstruction scheme is: 



1. 



Calculate d^' = ■^^2<:L^, 



2. Calculate " ^El ^'^''^ i ^'^-"l 

»V a, a, J 

for the sake of clarity, we haven 't considered the boundary situation, but we WilJ addrws this 
later. 

CgreHHT41 . Suppose biorthogonal filters (A, h, g, g] are from initial filters {h, h\ g\ g] by 
the lifting scheme (4.1a) or (4.2a). If the decomposiuon and reconstruction by filters 
\h. h\g\g\cznht accomplished only by integer calculation, such as Example 2. we also can. 
create a coircsporiding integer wavelet decomposition and reconstruction scheme which is very 
"close" to the original one by using niters {h, h, g, g] . Here the word "close" means that the 
difference of the two decomposition schemes is just some rounding error, and this rounding error 
will be corrected by, the integer reconstruction scheme. 

In fact, if {:;•'} and {d[ }are integer after (4.-3), we can calculate {ci } by 

c[ -cl'^J^J^dlA (4.7) 

instead of (4.4). Here Int(x) , as described in Secuon 2. is an arbitrary rounding up fiinction 
which satisfies x - 1 ^ Int(x) i x f l.It is obvious that (4.7) is very "close" to (4.4); and the exact 
reconstruction scheme can easily be obtained from 



(4.8) 



and (4.6). There will be a similar result, if the set of biorthogonal Filters [h. h. g. g) is obtained 
from the iniUal set of filters {h\ h. g. g°) by using (4.2b). 

We can now note, except for the example shown in the Lazy wavelet, (Example 2) most 
standard biorthogonal wavelet tran-sfonm cannot be performed directly by integer, even for one of 
the simplest wavelets, the Hoar wav<tei. However. If we properly choose the parameters «, and 
or, . and slightly change the uansform algorithms, such as Example 1 and Example 3, we c^n have 
a variaUon of the original biorthogonal wavelet transforms with re.spect to the set of filtere 
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[h, h°, g} (or {A^. g^}). On the other hand, the parameters {j, } should be also 
chosen carefully to guarantee that only addition and shift operations are needed by the algorithm. 
Another observation: if the set of filters {h, h, g,g] is obtained firom a set of filters 
g^} by the lifting scheme, and the set {A^ g. g^} is also obtained from a filter set 
{/»°. h^, g^, , we can repeatedly use Corollary 1 to get a "close" integer wavelet 
transformation. 



5. The Correction Method for Creating Integer Wavelet Transforms 

In this section, we will describe another approach for obtaining integer wavelets by using 
the so called Correction method. The motivation of thi.*; method is from the S+P transform, and 
we will now generalize this approach. Actually, the lifting scheme for generating biorthogonai 
wavelets can be considered as a special case of the correction method. From this method we can 
get some even complicated filters with fast decomposition and reconstruction algorithm. 

Suppose that we already have a simple integer wavelet transform, such as Examples 1 
through 3» the decomposition and reconstruction scheme of which can be formulated a.s follows; 

Reconsuuaion < = ' (5.2) 

Here. (5.1) and (5.2) can be the same as (4.3) and (4.6) or other algorithms. 

In general, after the above decomposition, one may not be satisfied with the result. There 
may still be some correlation among the highpass components because of the aliasing from the 
lowpass components, or the lowpass components do not cany enough of the expected infonnaUon 
from the original signal. Hence, wc could make an impyrovemcnt by putting some correction part 
on the highpass components or lowpass components. There are many ways to accomplish this. 
However, for the sake of the integer calculation, we prefer to use following conection method. 
For example, if we want to make a correction for the highpass part, the corresponding formula 
would be: 



d[ ^d['-lni(dc\)k^ -: 0, 1. 2,-\ (5.3) 
Here, dc\ is a correction quantity for d[ 

<ic\-^t, a,c\t + Z r^dt° • k - . • 0 , 1 . 2, • • (5.4) 
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and, {a, J.'^^ and ^- j^^ are given parameters which have been chosen for the user's purpose, 

such as reducing the redundancy among bighpass components or some other special requirement. 
We are not going to discuss how to choose these parameters, but one can refer to the references [3, 
5, 6} (or cluifMtion of thii proun. The only thing we need to mention ii, for (lie Mke of the 
integer calculation, any entries in both {a,- and should be rational numbers with 

denominators being powers of 2, 

From (5. 1), (3.3) and (5.4), it is easy to see the perfect reconstmction algorithm can be 

rfi-' = dj+lnl(dcjk=- -. m, m-I.m-2-". (5,3) 

combined with (5.2). 

As mentioned above, the Lifting scheme is a special condition of the Correction method. 
Examples 3 through 5 can also be considered as the examples of this method. We next give an 
example of the Correction nacthod which cannot be included in the group of Lifting scheme, and 
also which does not result in a closed form of compact .support for biorthogonal filters. 



Riramni- <t s^P tiansform [3], which is similar to using following a nalysis filters 



n 


-2 


-1 


0 


1 


2 


3 


^' 


0 


0 


1/2 


1/2 


0 


0 


z. 


n/16 


-1/16 


15/32 


-17/3i 


7/32 


-1/32 



While, the synthesis filters do not have compact support. However, the S+P transform can 
be implemented as follows: 
(a) Decomposition 

( 1 ) Take the decomposition step of Example 1 , that is, compute 



and 



/l.O 



=Int(-^)+CjV,. Jfc = 0,-..A^, -2, 



Int(— ~) + c^^ . if iV is an even number. 



if N is an odd number. 



(2) Correction Step: Define 5o - -I, Ji = 1 . 7' -1 and 
and now compute 
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</i-c/r-Int(lz£L); 

4 



(b) Reconstruction 
(1) Compute 



(2) If f/ is an even number, compute 



^1 



or, if A/ is an odd number, we have 



(3) Compute 



6. Boundary Conditions 

In the previous two sections, wc did not show how to get the integer, wavelet transform at 
the boundaries of signals. However, for all of the examples given above, the boundaries have 
been considered. TTiere are two issues dealing with boundary filtering if we use the L^ing scheme 
or the Ccrrection method to generate the integer wavelet transformations. The first is how to 
process the boundaries which occur in the surt-up wavelet iransfonnaiion.v The second is how to 
deal with the boundaries in the deductive formula If the boundaries in the startup wavelet 
trwisfoim have already been established, then those in the upgrading formula are easy to establish. 
In fact, for the Ufting scheme, the boundaries in both steps should be processed in the same way. 
While, for the Correction method, it is easy to see from (5.3)-(5.4) that one has more choices to 
process boundaries in the second step. Therefore, the only thing we need to discuss here is the 
process by which the boundaries in the start up wavelet transformations established. Assume 
we begin with compaa supported biorthogonal wavelets. 
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Suppose the original signal is {c^}^^. For creating integer biorthogonal wavelet 

transfonnations we can use the following symmetric extension [7]: 

(1) . If current biorthogonal filters have even length, we extend 
the boundaries of the signal as cf^ =^4^. k »l,2/s 

(2) . If the filters have odd length, we do the extension as c!, 2 cf , k » 1 ,2r • • 
Example 1 through 5 use the boundaries give above. In Example 6, the start up wavelet 

transforra uses the above boundarieji, but in the upgrading step, another boundary filtering is used. 
In addition, for arbitrarily sized images or signals, one can use the same technique which we 
described in the above examples to deal with this condition. 

7. Some Applications 

Before talking about any applications of the integer wavelet tnmsform given above, we fust 
prove that a nice pToperxy of precision preservation (PPP), which is similar to the one mentioned in 
Section 2, holds for both the Lifting and Correction upgrading technique. This property is very 
important for many applications. 

Lemma 7,1 Suppose that our integer wavelet transform starts with a pair of biorthogonal fillers 
with the PFP property discussed in Section 2, that is, (4.3) and (4.6) possess this property. Tbcn, 
the same property will be preserved in ihc whole-algorithm if we adopt the Uftii^j scheme to be 
the upgrading formula. 

In other words. Lemma 7.1 states if wc only use the working units with the same precision as the 
original signal or image to calculae the wavelet transform developed in Section 4. the equations 
(4.8) and (4.6) are still the backward operations of the cquation.s (4.3) and (4.7). 
Proof. Assume that we only use g bits to represent images or signals, say, the range of the pixel 
values is within [-2^' , 2^' - 1] . According to the hypothesis of the lemma, the equations (4.3) 
and its invert (4.6) have the PPP property. Therefore, what wc have to verify here is that the 
equation (4.7) and its inverse (4.8) can preserve the same property. We rewrite (4.7) and (4.8) as 
follow: 

ci=ci''+*,. (7.1) 

and its inverse 

c;' = c;-fr,. (7.2) 

Here, 
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In fact, the quanUty b, would have the same value in both (4.7) and (4.8) if we calculate it 
in the same way. On the other hand, if the working unit for b, is g biu, the machine will give b, 
another value, say ( - 2- S fr, < 2-^ ). where b, is not equal to in the sense of 
mathematics if the value of b, is beyond the interval [-2*-'. 2-' - 1] . However, b, will be the 
same in both (4.7) and (4.8). Therefore, the machine will automatically implement (7.1) and (7.2) 
as 

;/ -Td i ^ _i.o 



cr+i*. if .Z'-'Scl^^". <2«-', 

c:°-H^,-2*. ifcr + i,>2-'. (7.1m) 

^'^c^'+b, if c;%^,<-2'-'. 



and 

'.I 



<-b,. if .2*' <cl-^, <2*-', 

2*+ci-i,, if c\-b, <-2'-*. 

cl-b,-r, if c\^b,ir^. 



(7.2m) 



It is easy to see that (7.2m) is just the backward operation of (7. 1 m). which provides the evidence 
that the conclusion of this lemma is conecL 

It should be mentioned that the coefficients {cl ] obtained by (4.3) and a.Im) might not be 
the "real" wavelet coefficients using common sense. However, if we still use the working unit 
with q bits precision ar tbc rcconstixicUon step. (7.2m)*tnd (4.6) will give the exact original signal 
back. On the other hand, the coefHcienu {cj } still keep the most continuity of the "real" wavelet 
coefficients. Theiefore, when we repeat the decomposition step on {c| }. most small coefficients 
in its high frequency part {rf.' } will be almost the same a..; the "real" coefficients (within some 

rounding eiror). which aUows us to still take advantage of the "real" wavelet transform in image 
compression. 

A similar argument can show the same PPP property will hold for the integer wavelet 
transforms generated by the Correction mediod in Section 5. 

As we mcnUoned before, for many applications, the lossless image compression is as 
important as lossy compression. The integer wavelet transforms give the opportunity to compress 
without loss. It is also obvious that the integer wavelet algorithms can be used wherever ordinary 
wavelets are used, especially in signal and image compression. However, for most computers, the 
integer wavelet transform is much faster than the ordinary one and it uses much less mcmoiy. 
The following are some applications illustrating these types of transforms. 
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AppUcatioD 1. Lossless image compression 

As meationcd at the beginning of this paper, the integer wavelet transformation 
established by the techniques described in this paper can always be used for lossless image 
compression because of the reversible ability. Especially, wc can use the PPP property 
discussed in Lemma 7. 1 . We have used this wavelet lossless technology (WLT) for gray 
scale lossless image compression, and we have tried several images. For most natural 
images, the size of wavelet lossless compressed images is much smaller than corresponding 
GIF images. Figure 1 through 4 give some examples. Figure 1 is a standard image for 
compression. Figure 2 and 4 arc X-my images and Figure 3 is a two-value image but we tr«at 
it as a 8 bit gray scale image in order to compare with the GIF fomiat In fact, if we convert 
Figure 2 to a binary image, better result can be obtained by the ,IBIG technique. = 




(51U512) 
Fijarf 1. Comprcwion Raiio 
WtT: I.9:1/GIF: 1 05:1 




(512x512) 
Ficur* I. Compression Rstio: 

WLT4.5:l/GIF: 2.72:1 



"Visioneer may have come up with on 
against the paper blizzard. . .gets piles 
way to others throughout your compar 

Figure 3. Compression Ratio: WLT: 20.8/ GIF 17.8 (132x794) 
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Figure 4. Compression Ratio: WLT 3.8:1/GIF 1.98:1 (1232x1024) ^ 



Application 2. Large scale medical image compression 

Usually, 12 bits are used to rrpreseni one pixel in niedical images. In this situation, the 
values of the pixels vary from 0 to 4095. Such images require careful treatment when a 
transform coding method is used for compression. If we use ordinary biorthogonal wavelets, 
the range of the transform cocfTicicnts will expand to [-2'^ 2'*] when five levels of transform 
are used Therefore, a longer working unit has to be employed, which consumes significant 
computer resources. However, the integer wavelet technique developed in this paper will solve 
this problem. For example, if wc use the transforms given in Example 3. 5 and 6, the values 
of transform coefRcicnts will be linuted to the range of [-2'*. 2*^]. Even if wc do not use the 
PPP property for these wavelets. 16 bits for the working unit is sufficient for all computations. 

8* Conclusion 

This paper has shown the processes necessary in order to obtain a non-linear, integer, 
biorthogonal, or non-biorthogonal reversible wavelet n^sfomi suitable for signal or iniage 
processing. Wc have shown how such a transform can be obtained cither using tbe Lifting 
method, or the Correction method. For example, all interpolation wavelets can be modified to be 
corresponding integer wavelets without loosing any properties of original wavclew. In addition, 
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we have shown under certain conditions, the precision of the transform computation on the 
computer can remain at the same pircision of the data, thus reducing the need for additional 
computer memory during the transform computation. These arc extrtmely po^v«ful techniques 
when the target data are large images, or the requirements establish a need for speed. 

Although this paper establishes the structure for the integer tmisform based upon the 
biorthogonal wavelet or some non-bioirhogonal wavelet, we do not imply the examples in this 
paper are necessarily the best wavelets for any particular applicition. However, we do claim if 
one is going to use such a technique, the ideas suggested in this paper wiU provide the best 
implcracntation. 
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Compressor V3.0 




Armouadsg, Lightning Strike™ Image Compressor (LSIQ version 3,0, a Windowi 95 tool that 
compresses still images from 50:1 to 200:1 using INFINlTRON'i proprieUry wavelet technology. LSIC 
is a versatile, tMsy to use tool for Web and Graphic designers that can handle a wide variety' of digital 
image formats, and it fftdndes filters tad convenient web tools. Images can be compressed 3 to 5 times 
more than JPEG, while maiAtaioing similar or better image fidelity. Images can be viewed in 2 to 4 
seconds over the Inurnet rather than 10 to 20 seconds for images compressed under JPEG. This has 
enormous benefits for reducing bottlenecks on coq>orate networks and the web* snd in addition, requires 
less ftorage space. 




liglitnliB Strlli Fcituret 



Compressiom. Images can be compressed as high as 
200:1 using wavelet lecbnulogy. 

Comprcssioa Cootrol. An EASY mode allows the user 
tu compress images with minimal input, requiring only a 
deci^OQ between wor% quality or more compresifiixi. An 
ADVANCED mode enables ihc user lo selea; 1) image 
fils yi>«. 2) compressioD ratio. 3) PSNR. or 4) master 
Icvd. Web designers will like the one step pnxess to 
control the size (}f their image files, thus insuring the 
j^eod an image may be viewed uc a brow 

Non Uoiforoi Compretsloi. Re^ons of an imAge can 
be selected for leas COTiprcssion to preserve a highff 
image quality while the rest of the image is comivessed 
10 the specified comprcjwion ratio. In this way importaDt 
parts of a picture maintain crucial details while thopver^^ 
all picture file can be made aa small as possible - 



Post Rcconstructiou Filters. FtUerii arc available to 
enhance the rcconjuructed image. At compre.'Viioo time 
the user can preset a control to have these f^Ucrs operate 
automatically durinj^ rcconstructiotL The filters include; 
quality imprtsvemeot, sharpen (edge enhancement), 
smcx^ing, and brighten. 

Traosparcadcs. Tbe user vdll hove the ability to set 
pixels txaosparent » that a color in the background 
(already on the page) can be seen through the picture. 
ThLi is useful for creative web site developers. 

Progressive Compreastoa. An image can be compressed 
so that when it is viewed it will appear quicWy, first with 
low resolution, and then progressively building up in 
detail as ii in downloaded. This insure.** the viewer does 
not lose iniae.li while tbe^ unage is duwnloaded ^ 
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Cow, Bitxsap, No Coopres&ioo 




Cow, LigbtniBC Stnke, &5:1 compressioa 



abfuiiiFiiiTiii. lie. 

FouTxkd in 1992, INFINTTRON, Inc. is a private 
cofnpany ihas apecializo in the design and marketiD;; of 
high quality image and video compresairn solutiooa for a 
wide anay of markcw. INFINITRON i« based in 
Vancouver, BC with labs in Rcgina. Saskatchewan and 
Dcntoo Texas. 




Cow, JPEG 1^5:1 CconpressioQ 



Apidcadm Uittniai Srilc iMigt Maprtniei 

images on the Web 
Pb«o Stock 
Data Warehousing 
Catalogues 
Video Games 

CDs (Hncyclopcdias, Museums. Science, and Medicine) 
Archives (An, History, Genealogy) 
Medical Imaging 

Perlirninci 

Encode lime typically less than 3 seconds for a 320 X lAQ ^ 
pixel, 24 bit color image on a 133 MHz. Pentium with 16 MR 
RAM- 

■inimuB Rccimmeidtd tisttn 

Windows 95/NT OS 

Pentium 100 MHz, 8 MR RAM 

2 MB fur program files 

10 MB plus to swap image files 



JtacBlwy linimilN Pniicti 

Na»:apc Navigator Plug-In 
Java Applet 
ActiveX Control 
Web Site Image Convener 
Lifihtoing Surike SDK 
GML Banner GenftrAior 




iX)wnload a FREE deno version of Lightning Scikc Windows Cwnpresw from: 

www infmitroocom 

INHNTTRON 3401 East University, #104, Dcntoo, TX, 76208 
USA Office Tel: 817.484.1 165 FAX: 817484.0588 

"INFINITRON ao*" Flr-^ 1 199 W. Hastings, Vancouver,3C^V6K-3T5 
Canada Office Tel: 604.688.9789 FAX: 604.688.9798 ^^s^ 
March 1 W7 Crmvnchl INFTNTTRON Roiearch Inwmaticm*!. Tnc- All Riehu Rc»crvorf 
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Announnng, Ughtning Strike Imm^t Compreswr (LSIC) version 3.0, a Windows 95 tool that 
compresses rtfll images from 50:1 to 200:1 using InCaitron's proprietary Wavelet technology. LSIC is a 
versatile, eagy to use tool for Web and Graphic designers that can handle a ^de variety of digital image 
formats and includes filters, and convenient web tools. Images can be compressed to Mes 4 times smaller 
than JPEG, while mainuining similar or better image fidelity. Images can be viewed in 2 to 5 seconds 
over the internet rather than 10 to 20 seconds for images compressed under JPEG. This has enormous 
benefits for reducing bottlenecks on corporate networks and the web, and in addition, requires less 
storage space. 




llBHtiini Strike Features 



CompresiioB. Images can be compressed to'uu«r 
usin^; Wavda Technology. o-^ ' i 

Compres&ioB Control. An EASY mode allows the user 
to comprcjw images with minimaJ input, requiring only a 
decision between more quality or more compression. An 
ADVANCED mode enable? the usa to select; 1) image 
flic size, 2) compression ratio. 3) PSNR, or 4) master 
level. Web disigncrs will like the one step process to 
control ilie size of their image files, thu.'i injuring the 
?!peed an image m^y be viewal on a browsa. 

Non Uuifornj Compression. Regions of an image can 
be itelccied for less compression to preserve a hiiter 
imA^ quality while ihc xcsx of the image i.« comprejwcd 
to the specified eornpresijiun ratia Jn this way impfiriacl 
pans of a picture miinuiin crucial details while , the over ■ 
all picture fi J c can be made as small as possible 



PoU Reconstnictioo Filters. Filter.'! are available to 
cniviDCB the reconstructed image. At compression time 
the user can preset a control to have these filters operate 
Quiomaiically during reconstruction. The filters include; 
quality improvement, sharpen (edge enhancement), 
Sfmoothing. and brighten. 

Transpareades. The user will have the ability to set 
pixels transparent so that a colur in the backi;round 
(already on the page) can be seen through, the picture. 
This is u.<5cful fowr creative web site developers. 

Progressive Compression. An image can be compressed 
50 that when it is viewed it will appear quickly, first with 
low resolution* and then progressively building up in 
detail as it is downloaded. This insures the viewa does 
not loose iuTcresi- whilcnhe image is download ed."'^"* 
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Lcnna .bmp. No Compression 




Ahout L^n^^aTRON. inc. 
Fo'^dcd in 1992, INFTN7TR0N, inc. is a private 
company ihai spccialim in the design and maricaing of 
lii£h quality image and video compression soluiions fnr a 
wide array of markcw. INFINITRON is baswl in 
Vancouver, BC with labs in Regina, Saskatchewan and 
DciOD Texas. 




Lcnna. JPEG 115:1 Compression 



A»liflc«(lMs UgUilna nrfki taige CinvrtssiMi 

Images on ihc Web 
Photo Stock, 
Data Warehousing 
Caialoqucs. 
Video Games 

CDs (Encyclopedias, Museums, Science, and Mediuinc) 
Archives (Art, History, Geneoly^) 
Medical Imujpny 



Pirrirmiice \ 

Encode time typically Icsii than 3 seconds for a 320 X 240 
pixtJ. lA bit color image oo a 133 MHt:. Pentium with 16 MH 
RAM. 

Mlnlnun Recanmiidti Syitin 

Windowi? 95/NT OS 

PcDtium 100 MH7^ 8 \4B RAM 

2 MB for progTiim files 

10 MB plus to swap image files 

tacliinrlNHlimoiPnducu 

• Netscape Navigator Plug-in 

• Java Applet 

• ActiveX Conu'ol 

• Web Site Image Coovcnff 

• Lightning Strike SDK 

• HMT. Runner Hoieraior 




Download a FREE demo veriion of Lighming Suikc Windows Compressor from: 

www.infinitron.com 

INFTNTTRON 3401 East University, #104, Denton, TX, 76208 
USA Office TcJ: Kl 7.484.11 65 FAX: 817.484.0588 

JNFINTTRON 1 C^-^Klf^l ] 99 W. Hastings. Vaocouve. HCrV6E-3T5 - • : . • V 

Cinada Office Tci: 604. 9789 FAX: 604.688.9798 
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The Lightning Stnke CompnsMor is a Windows tool Ihit comprwe* still images from ■ wide 
variety ofdigital image fomuts using Infinitron's proprietary Wavelet algorithm. Images can 
be compressed to fiUa 5 Umcs smaller than J?EC. while mainuining similar or better image 
fidelity. Images can be viewed in 1 or 2 .vxonds over the internet rather than 10 to 20 
seconds for images comppwsed under JPFG. This has enormous benefits for transmitting 
over corporate networlcs or the web, and in addition, saves space required for storing all those 
imagci. 

Lightning Strike is i collection of toots in a user friendly environment Two levels of user 
control are offered, one quick and easy for most applications, the other a master level for the ' 
advanced user who wishes to control parameters to maximize image quality. 

The compression approach used by Lightning Strike is based upon integer wavelets. This 
technology is acknowledged by leading experts » a superior comprc»ion technique u 
conipaiied lo discrete cosine minsfonn used in JPEG. 



ireS:.:-.:- 

Image Compression Options and Control 

Compression Technique Options 

Both Infuiitfon's Wavelet Compression and other frequently used compiwion metJwds are 
included in the product so users need only have Lightning Strike on their work station to 
^orm bI] image compressions. Images can be compressed to Wavelet, JPEG. PNC, and 

Compression Quality Versus Speed Options 

The user can select one of two encoding processes that trade quality for speed of compression 
and eaae of uae. With the -Advanced- option selected, the optimum compitwion parmmetcm 
w »et^by the user to give tbc best possible images for selected comprwsion r«io. With 
^^^ected you get the £>stest compreuion without having to know details of parameter 

Comprwsion Ratio Control 

Jecompressed^^^^ file lire or compression ratio may be specified rather dian the quality 
ftctor. This enables a web designer to control the size of their image files or the speed an 
image may be viewed, in a one step process. 

Region of Interest Focusing 

Rcgioru of an image can be selected for lc» compn^sion to preserve a higher image quality 
while the rest of the image is compressed to the specified compression ratio. In this way 
important parts of a picture maintain cnjcial deuils while the over all picwre file can be made 
as small as possible. This is also known as Non-Uniform Compression. 
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Sptft and Merge 

Very iirjc images, which could not otherwise be compressed due lo their Urge size, can be 
split into smaJlcr images and compressed individually. Hiis process has the side advantage of 
using RAM more eflteiveiy speeding time for comprcsxion. The «plii images can be 
reassembled using the merge aspect of tlie feature. 

Post Reconstruction Fiilers 

Filten arc available to enhance the reconstrucicd image. At compression time the user can 
prwat a control to have these filtCT operate automaiically during reconstruction. The filters \ 
include; quality improvement, shwpen (edge enhancement), smoothing, and brighten. 

Transparencies 

The user will have the ability to set pixels transparent so that a color in the background 
(already o;i the page) can be seen ihrougj) the picture. This is useful for creative web site 
dcvelopcnt. This givc« the ability lo display pictures other than the rectangular shape allotted 
on ihe web page, i.e. circles, polygons etc.. Also, designers often use this feature for 
shadowing, letters and objects. 

Progressive DecompfBssion 

An image can be compressed so that when it is viewed it will appear quickly, first with low 
resolution, and then progressively building up in detail as it is downloaded. This insures the 
viewer doc» not loose interest while the image is downloaded. 
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Picture of Una, n^ith no Compression (512X512 Image) 
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PNG FUo Structure 

The cotnprwMd image* trc «ored in file formni compliant with the PNC lUndaid. In the 
future this flic formar will replace tUe GIF format us«d today. 

Batch Comprtss 

The u»«r is able to comprcM many images at once by adding or deleting image filejt (or path*) 
to a list box. 

Imag« Statistics 

The compressor stores Image sutistics on each compressed image which may be viewed by 
the user. The following infomtaiion is provided; image dimensions, compression ratio, file 
sizes, MSE, PSNR, maximum pixel difference, compression and decompression limes. 




Encoding and Dacoda Tima. 

The lypicai tinvi to encode or decode a 320X240, 24 bit color image is I second on a Pentium 
runnhtg at 135 NSIz with 16 Meg RAM. 

Minimum Racommandad Systam 

The minimum system irquircmcnts for an IBM PC Compatible are: 

Hard Disc Drive 2 Mbytes free for program files, 10 Mb plus to swap 
image files. 

Operating System MS Windows 3. 1 (Wiii32y 95/ NT 
RAM 8 Mbytes 

This software U also available on the Apple MAC, Solaris, and UNIX platforms. 
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Andllaiylnfinlfinon Products 
Netocapa Navigator Plug-in 

Netscape plug-ins are available for the Mac68k, Mac ?PC. Widows 3.1 and 95/NT. 
Java Applot 

Java Appl€t5 arc available for the Mac6« k. Mac PPC, Widows 3. 1 and 95/NT. 
ActiveX CoEitrol 

The Lightning Strike decompression sofiwarr h available for such applications as Microsoft's 
Internet E^tplortr, as an ActiveX control. 

W»b Sitt ImagQ Converter 

This utility will automate the conveaion of web pages from JPEG to the Lightning Strike 
fonnat The utility searches an HTML file and replicates it replacing any JPEG image tap 
with Lightning Strike tags and converting the JPEG image files to Lightning Strike. The 
utility can follow link tagi to recursively convert and replicate an entire web site or sub- 
section of a web site to the Lightning Strike format This utility will be available for 
Windows NT and most flavors of the UNIX operating system. 



Lightning Strike Software Developers K\X 

Using the SDK, a developer ean integrate the highly efficient Lightning Strike module 
libraries into their own applications. 



Download a FREE demo version of Lightning Strike Windowi 
Compressor from: www.infinitron.com 
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Windows Image 
Compressor V3.0 



Announcins, Lightning Strik*~ Image Compressor (LSIC) version 3.0. . Windows 95 tool that 
compress^ stJI images from 20:1 to 200:1 using INHNITRON's proprietary wavelet technoloey LSIC 
.s a versatile, easy to use tool for Web and Graphic designers that can handle a wide variety off gital 
.mage formats^ and it includes filters and convenient web tools. Images can be compressed at ratios well 
in excess of pr«ent JPEG raUos while maintaining comparable image fidelity. This translates to much 
shorter image down load time on the web. This has enormous benefits for reducing bottlenecks on 
corporate networks and the web, and in addition, requires less storage space. 
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CompressioQ. Uses a proprietary integer wavelet. 

Comprcsjion Control. An EASY mode allows the user 
to compress images with minimal input, requiring onJy a 
decision between quality and compression. An 
ADVANCED mode enables the user to select; 1) image 
file size, 2) compression ratio. 3) PSNR, or 4) master 
level for professionals where every parameter can be 
altered. We also provide the highest, wavelet lossless 
compression for users wishing this capability. Web 
designers will like the one step process to control the size 
of their image files allowing control over the delivery 
time of an image over a network. 

Noa Uniform Compression. Regions can be seleaed for 
,l^.:S95^P^-*o" ^0 preserve image quality while the rest 
compressed to the specified compression 
ratio.. In this way. important parts of a picture maintain 
cnicial details while the over ail picture file can be made 



as small as possible 

Port Recoortruction Filters. Filters are available to 
enhance the reconstructed image. At compression time 
the user can preset a control to have these filters operate 
automatically during reconstruction. The filters include; 
visual quality improvement, sharpen (edge 
enhancement), smoothing, and brighten. 

Transparencies. The user will have the ability to set 
pixels transparent, so that a color in the background 
(already on the page) can be seen through the picture. 
This is useful for creative web site developers. 

Progressive Compression. An image can be compressed 
so that when-'ii is viewed it will appear quickly/ first with 
low resolution, and then progressively building-up in 
detairas it is'downloaded. This insures the viewer does 
not lose interest while the image is downloaded. 



I N F I N ITR O N 




Cow. Lightning Strike, 85:1 compression 



About IMFINITieiJic. 

Founded in 1992, INFINrrRON, Inc. is a private 
company that speciiiizes in the design and marketing of 
high quality image and video compression solutions for a 
wide array of maricets. INFINITRON is based in 
Vancouver, BC with offices in Regina, Saskatchewan and 
Demon Texas. 



]Vjndoivs Image 
Compressor V3,0 




Cow, JPEG 85:1 Compression 



AfspflctdtDS !ir Uattnlag Strtta Isuoi CMumsfUo 

• Images on the Web 

• Photo Stock 

• Data Warehousing 

• Catalogues 

• Video Games 

• CDs (Encyclopedias, Museums, Science, and Medicine) 

• Archives (An, History, Gcncaloa^:) 

• Medical Imaging . X 

Pirtiraaoci 

Encode time tvpically less than 2.5 seconds for a 640 X 480 
pixel, 24 bit color image on a 133 MHz. Pentium with 16 MB 
RAM. Decode time is less than .75 seconds. 

HlalDUB lacinaiaBdai Systaa 

Windows 95/NT OS 

Pentium 100 MHz, 8 MB RAM 

2 MB for program files 

10 MB plus to swap image files 

tacuanr mHNmoi prt docu 

• Netscape Navigator Plug-in 

• Java Applet 

• ActiveX Control 

• Web Site Image Convener 

• Lightning Strike SDK 

• GML Banner Animation/Compression 

• Black and White Image Compression 




Download a FREE demo version of Lightning Strike Windows Compressor from: 

ww\v'.infi nitron.com 

INFINITRON 3401 East University, #104, Denton, TX 76208 
:USA0f5cc Tel: 817.484.1165 E^AX; 817.484.0586 

INFINITRON 10* Rr 1 19*9. W. Hastings. Vancouver, BC, V6E- 3T5™ 
Canada Office Tel: 604.688.9789 FAX; 604.688.9798 



This Page is Inserted by IFW Indexing and Scanning 
Operations and is not part of the Official Record 



Defects in the images include but are not limited to the items checked: 

□ BLACK BORDERS ' ' 

□ IMAGE CUT OFF AT TOP, BOTTOM OR SmES 

□ FADED TEXT ORDRAWING 



(^V^<:<;OLOR OR BLACK AND WHITE PHOTOGRAPHS 

□ GRAY SCALE DOCUMENTS 

□ LINES OR MARKS ON ORIGINAL DOCPMENT 

□ REFERENCE(S) OR EXHrorr(S) SUBMITTED ARE POOR QUAtriY 

□ OTHER; " • ■ 

IMAGES ARE BEST AVAILABLE COPY. 
As rescanning these documents will not correct the image 
problems checked, please do not report these problems to 
the IFW Image Problem Mailbox. 



BEST AVAILABLE IMAGES 



Defective images within this document are accurate representations of the original 
docimients submitted by the applicant. 



□ BLURRED OR ILLEGIBLE TEXT OR DRAWING 



□ SKEWED/SLANTED IMAGES 




